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^  ABSTRACT 

The  line  array  solution  of  Part  I  is  used  in  convolution  with  the 
primary  pattern  to  calculate  the  field  of  collimated  and  spherically 
diverging  sources  of  arbitrary  transverse  distribution. 
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INTRODUCTION 

In  Part  I  (reference  (1))  the  contour  integration  method  was  used 
to  calculate  the  field  of  a  line  array.  On  the  axis  the  solution  has  a 
logarithmic  singularity.  Finite  aperture  axial  levels  were  calculated 
for  two  special  cases:  a  cylindrical  and  a  conical  source  in  which  the 
transverse  densities  are  constant  within  the  aperture  and  zero  outside. 
A  "complete"  field  solution  for  each  case  is  now  obtained  in  which  the 
transverse  distributions  can  be  arbitrary.  The  solution  is  in  the  form 
of  a  convolution  between  the  line  array  pattern  and  the  source  pattern. 
This  is  equivalent  to  summing  the  fields  of  a  bundle  of  line  arrays  of 
appropriate  weights. 


*(R) 


CYLINDRICAL  SOURCE 

The  volume  integral  for  the  cylindrical  source  can  be  written 
Q0  exp(-aL-ikL) 
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Where  t|>  is  the  velocity  potential,  Q0  is  the  linear 
Is  the  range  vector  to  the  field  point  and  aj2  =  a2 
where  V'is  the  cylindrical  angle  of  the  radius  a', 
normalized  transverse  distribution  function;  i.e. 


a ’da'  D(a 1 ,y 1 ) 
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source  density,  R 
+  a'2  -  2aa'  cos(y'-y) 
D(a' ,y ' )  is  the 
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Let  Y  *  L  i|i(L,a)exp(ikL)/S0  where  S0  =  Q0/4to  is  the  Westervelt  source 
strength.  Equation  (1)  can  then  be  written  in  dimensionless  form  as 
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D(a 1 »y ' )  Y(u0,v0) 


From  reference  (1) 


Y(u0,v0) 
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dz  exp(-z) 
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Is  the  line  array  solution  and  where  z  +  -u0  +  ivQ  *  -2oL  +  ikai2/2L 
and  B2  «  41u0v0  for  L»a. 


Equation  (4)  Is  easily  evaluated  along  the  contour  z  *  t  -  u«  +  1vQ. 
Equation  (3)  Is  then  the  two  dimensional  convolution  of  Equation  *4)  ana 
OU'.y'). 


DIVERGING  SOURCE 

In  the  spherical  case  let  Y  *  R^ (R)exp(i kR)/S0 .  Equation  (3) 
becomes 


Y  = 


2ir  it 


d*'sin*'  D(*\y')  exp(iv0)  Y(u0,v0) 


where  D  is  the  normalized  density  function; 
2  IT  IT 


d*'  sin*'  D(*',y')  5  1 


i.e. 
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and  where  Y(u0,v0)  is  again  given  by  Equation  (4).  For  small  angles 
u0  «  2aR,  v0  *  uo0j2  and  0i2  =  e2  +  9'2  -  296  'cos(Y  '-y)  where 
0'  *  * '/*n i  *0  being  the  characteristic  Westervelt  angle  of  reference 
(1).  At  long  ranges  Y(u0»v0)  approaches  the  Westervelt  pattern  except 
for  the  singularity  at  =  0.  Some  results  of  convolution  in  this  limit¬ 
ing  case  are  reported  in  reference  (2). 

Equation  (5)  was  programmed  for  numerical  evaluation  using  the 
conical  beam  approximation.  The  singularity  at  0j  =  0  was  avoided  by 
using  small  finite  initial  values.  The  results  for  various  values  of 
*uq  ■  2oR  are  shown  in  Figures  1  -  3.  The  scaled  conical  angle  is 
given  by  e0  =  *0'/*p  where  *0'  is  the  half  width  of  the  conical  beam. 

The  abscissas  are  the  relative  angle  9  *  */*0  as  in  reference  (1).  The 
ordinates  are  20  log^g  |Y|.  The  dotted  curve  is  the  Westervelt  pattern. 

For  90  *  0  the  conical  beam  is  a  delta  function  and  the  pattern  has 
a  logarithmic  singularity  given  by  (u0«l) 


Y(u0,v0)  ♦  -uQ  exp(-u0)  In  |v0| 

0  +  0  (7) 
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For  90  t  0.  Y  rapidly  approaches  the  finite  axial  value  for  e  <  e0. 
For  e0  >  l  the  curves  cross  the  e0  ■  0  curve  and  approach  It  from 
above.  For  e0-*»  the  e0  ■  0  curve  acts  like  a  delta  function  with 
respect  to  the  conical  beam.  It  Is  clear  that  only  for  e0»l  does  the 
actual  shape  of  the  source  pattern  have  any  significant  effect  on  the 
result. 

Figure  4  compares  the  results  with  experiment  where  e0  =  0.4.  For 
2aL  *  1  the  agreement  Is  good.  The  -3dB  beamwldth  Is  only  %  the 
Westervelt  beamwidth.  For  2aL  =  0.25  the  experimental  axial  value  is 
again  roughly  3dB  low  as  It  was  In  the  calculations  of  reference  (1). 
Misalignment  of  source  and  receiver  is  one  possible  explanation  of 
this  error. 


SATURATION 

The  effects  of  saturation  may  be  included  as  in  Part  II  (reference 
(3))  by  multiplying  the  integrand  of  Equation  (4)  by  the  taper  function 
T2(z).  The  solution  for  the  cylindrical  and  spherical  cases  may  also 
be  combined  to  give  a  reasonably  good  approximation  for  all  ranges.  At 
moderate  range  the  spherical  part  alone  should  give  a  sufficiently  good 
approximation  if  the  result  is  multiplied  by  the  source  aperture 
pattern  (reference  (4)). 
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Figure  3.  Beam  Pattern 


Figure  4.  Comparison  of  experimental  and  predicted  beam  patterns. 


